This Rmarkdown script (and corresponding TAB-separated CSV input data file InfoRateData.csv and the resulting HTML document) contain the full analysis and plotting code accompanying the paper Human languages share an optimal information transmission rate.
For more information on the data, please see Oh (2015). There are in total 17 languages (see the Table below).
The oral corpus is based on a subset of the Multext (Multilingual Text Tools and Corpora) parallel corpus (Campione & Véronis, 1998) in British English, German, and Italian. The material consists of 15 short texts of 3-5 semantically connected sentences carefully translated by a native speaker in each language.
For the other 14 languages, two of the authors supervised the translation and recording of new datasets. All participants were native speakers of the target language, with a focus on a specific variety of the language when possible – e.g. Mandarin spoken in Beijing, Serbian in Belgrade and Korean in Seoul. No strict control on age or on the speakers’ social diversity was performed, but speakers were mainly students or members of academic institutions. Speakers were asked to read three times (first silently and then loudly twice) each text. The texts were presented one by one on the screen in random order, in a self-paced reading paradigm. This way, speakers familiarized themselves with the text and reduce their reading errors. The second loud recording was analyzed in this study.
Text datasets were acquired from various sources as illustrated in the Table below. After an initial data curation, each dataset was phonetically transcribed and automatically syllabified by a rule-based program written by one of the authors, except in the following cases:
Additionally, no syllabification was required for Sino-Tibetan languages (Cantonese and Mandarin Chinese) since one ideogram corresponds to one syllable.
| Language | Family | ISO 639-3 | Corpus |
|---|---|---|---|
| Basque | Basque | EUS | E-Hitz (Perea et al., 2006) |
| British English | Indo-European | ENG | WebCelex (MPI for Psycholinguistics) |
| Cantonese | Sino-Tibetan | YUE | A linguistic corpus of mid-20th century Hong Kong Cantonese |
| Catalan | Indo-European | CAT | Frequency dictionary (Zséder et al., 2012) |
| Finnish | Uralic | FIN | Finnish Parole Corpus |
| French | Indo-European | FRA | Lexique 3.80 (New et al., 2001) |
| German | Indo-European | DEU | WebCelex (MPI for Psycholinguistics) |
| Hungarian | Uralic | HUN | Hungarian National Corpus (Váradi, 2002) |
| Italian | Indo-European | ITA | The Corpus PAISÀ (Lyding et al., 2014) |
| Japanese | Japanese | JPN | Japanese Internet Corpus (Sharoff, 2006) |
| Korean | Korean | KOR | Leipzig Corpora Collection (LCC) |
| Mandarin Chinese | Sino-Tibetan | CMN | Chinese Internet Corpus (Sharoff, 2006) |
| Serbian | Indo-European | SRP | Frequency dictionary (Zséder et al., 2012) |
| Spanish | Indo-European | SPA | Frequency dictionary (Zséder et al., 2012) |
| Thai | Tai-Kadai | THA | Thai National Corpus (TNC) |
| Turkish | Turkic | TUR | Leipzig Corpora Collection (LCC) |
| Vietnamese | Austroasiatic | VIE | VNSpeechCorpus (Le et al., 2004) |
The data is structured as follows:
| Lng | # spkrs | % fem | # age | mean(age) | sd(age) | actual ages |
|---|---|---|---|---|---|---|
| CAT | 10 | 50 | 10 | 35.4 | 9.2 | (21, 28, 28, 29, 31, 39, 42, 42, 44, 50) |
| CMN | 10 | 50 | 9 | 23.1 | 4.5 | (19, 19, 19, 19, 24, 24, 25, 28, 31) |
| DEU | 10 | 50 | 0 | NaN | NaN | () |
| ENG | 10 | 50 | 0 | NaN | NaN | () |
| EUS | 10 | 50 | 10 | 28.0 | 4.9 | (19, 22, 26, 27, 28, 29, 30, 31, 32, 36) |
| FIN | 10 | 50 | 10 | 33.2 | 11.0 | (16, 22, 26, 28, 30, 35, 37, 41, 45, 52) |
| FRA | 10 | 50 | 10 | 32.5 | 7.7 | (24, 25, 25, 27, 28, 36, 36, 37, 41, 46) |
| HUN | 10 | 50 | 10 | 39.3 | 15.8 | (17, 27, 27, 31, 33, 39, 42, 51, 57, 69) |
| ITA | 10 | 50 | 0 | NaN | NaN | () |
| JPN | 10 | 50 | 10 | 30.6 | 12.8 | (20, 20, 21, 22, 22, 28, 29, 40, 51, 53) |
| KOR | 10 | 50 | 10 | 28.6 | 10.6 | (16, 19, 19, 19, 28, 31, 33, 35, 36, 50) |
| SPA | 10 | 50 | 10 | 33.7 | 10.1 | (21, 22, 26, 28, 30, 32, 42, 42, 44, 50) |
| SRP | 10 | 50 | 10 | 30.6 | 7.8 | (19, 21, 23, 30, 31, 32, 34, 34, 38, 44) |
| THA | 10 | 50 | 10 | 30.1 | 5.7 | (23, 23, 27, 28, 30, 31, 31, 32, 33, 43) |
| TUR | 10 | 50 | 7 | 32.6 | 7.2 | (24, 25, 30, 31, 37, 37, 44) |
| VIE | 10 | 50 | 6 | 27.2 | 4.1 | (21, 25, 26, 28, 31, 32) |
| YUE | 10 | 50 | 10 | 22.0 | 1.5 | (20, 20, 21, 21, 22, 22, 23, 23, 24, 24) |
NS: exploratory plots.
mean=101.196, median=100, sd=24.703, CV=0.244, min=49, max=162, kurtosis=2.484, skewness=0.177.
SR: exploratory plots.
SR per speaker.
SR by Sex and Age across Languages.
SR by Sex, Age and Language.
SR by language.
mean=6.631, median=6.777, sd=1.148, CV=0.173, min=3.589, max=9.492, kurtosis=2.408, skewness=-0.168.
ShE and ID: exploratory plots.
ShE vs ID.
Pearson's product-moment correlation
data: tmp1$ShE and tmp1$ID
t = 2.0326, df = 15, p-value = 0.06019
alternative hypothesis: true correlation is not equal to 0
95 percent confidence interval:
-0.02052914 0.77274779
sample estimates:
cor
0.4647009
Spearman's rank correlation rho
data: tmp1$ShE and tmp1$ID
S = 451.88, p-value = 0.07259
alternative hypothesis: true rho is not equal to 0
sample estimates:
rho
0.4462208
Paired t-test
data: tmp1$ShE and tmp1$ID
t = 11.635, df = 16, p-value = 3.213e-09
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
2.158040 3.119607
sample estimates:
mean of the differences
2.638824
ShE:
mean=8.621, median=8.69, sd=0.904, CV=0.105, min=6.07, max=9.83, kurtosis=4.665, skewness=-1.122.
ID:
mean=6.009, median=5.56, sd=0.883, CV=0.147, min=4.83, max=8.02, kurtosis=2.53, skewness=0.747.
ShIR: exploratory plots.
ShIR per speaker.
ShIR by Sex and Age across Languages.
ShIR by Sex, Age and Language.
ShIR by language.
IR: exploratory plots.
IR per speaker.
IR by Sex and Age across Languages.
IR by Sex, Age and Language.
IR by language.
ShIR:
mean=56.709, median=57.207, sd=9.35, CV=0.165, min=32.772, max=89.235, kurtosis=2.444, skewness=0.079.
IR:
mean=39.153, median=39.13, sd=5.097, CV=0.13, min=25.631, max=60.692, kurtosis=3.622, skewness=0.325.
SR vs ID
Pearson's product-moment correlation
data: info.rate.data$SR and info.rate.data$ID
t = -45.329, df = 2286, p-value < 2.2e-16
alternative hypothesis: true correlation is not equal to 0
95 percent confidence interval:
-0.7090031 -0.6658066
sample estimates:
cor
-0.6880138
Spearman's rank correlation rho
data: info.rate.data$SR and info.rate.data$ID
S = 3393600000, p-value < 2.2e-16
alternative hypothesis: true rho is not equal to 0
sample estimates:
rho
-0.6999614
| Level-1 factor (f) | ICC |
|---|---|
| Text | 0.30 |
| Language | 0.54 |
| Speaker | 0.00 |
| Level-1 factor (f) | ICC |
|---|---|
| Text | 0.01 |
| Language | 0.61 |
| Speaker | 0.30 |
| Level-1 factor (f) | ICC |
|---|---|
| Text | 0.01 |
| Language | 0.55 |
| Speaker | 0.35 |
| Level-1 factor (f) | ICC |
|---|---|
| Text | 0.02 |
| Language | 0.35 |
| Speaker | 0.50 |
| model | AIC | BIC |
|---|---|---|
| 1 + (1 | Text) + (1 | Language) + (1 | Speaker) | 2128.05 | 2156.73 |
| 1 + (1 | Language) + (1 | Speaker) | 2398.52 | 2421.46 |
| 1 + (1 | Text) + (1 | Speaker) | 2257.73 | 2280.68 |
| 1 + (1 | Text) + (1 | Language) | 4717.5 | 4740.45 |
| 1 + Sex + (1 | Text) + (1 | Language) + (1 | Speaker) | 2120.19 | 2154.6 |
| 1 + Sex + (1 | Language) + (1 | Speaker) | 2390.45 | 2419.13 |
| 1 + Sex + (1 | Text) + (1 | Speaker) | 2257.24 | 2285.92 |
| 1 + Sex + (1 | Text) + (1 | Language) | 4583.16 | 4611.83 |
We consider here the full model SR ~ 1 + Sex + (1|Text) + (1|Language) + (1|Speaker).
| model | AIC | BIC |
|---|---|---|
| 1 + (1 | Text) + (1 | Language) + (1 | Speaker) | 10255.96 | 10284.63 |
| 1 + (1 | Language) + (1 | Speaker) | 10522.63 | 10545.57 |
| 1 + (1 | Text) + (1 | Speaker) | 10305.76 | 10328.71 |
| 1 + (1 | Text) + (1 | Language) | 12886.67 | 12909.62 |
| 1 + Sex + (1 | Text) + (1 | Language) + (1 | Speaker) | 10244.07 | 10278.48 |
| 1 + Sex + (1 | Language) + (1 | Speaker) | 10510.53 | 10539.21 |
| 1 + Sex + (1 | Text) + (1 | Speaker) | 10298.66 | 10327.34 |
| 1 + Sex + (1 | Text) + (1 | Language) | 12743.83 | 12772.51 |
We consider here the full model IR ~ 1 + Sex + (1|Text) + (1|Language) + (1|Speaker).
We will use a Gaussian distribution (with fixed or modelled variance).
******************************************************************
Summary of the Quantile Residuals
mean = -7.129342e-05
variance = 1.000437
coef. of skewness = 0.04303611
coef. of kurtosis = 3.710291
Filliben correlation coefficient = 0.9979187
******************************************************************
Deviance= 1177.818
AIC= 1576.945
******************************************************************
Summary of the Quantile Residuals
mean = 0.002001853
variance = 1.000433
coef. of skewness = 0.03014488
coef. of kurtosis = 2.864355
Filliben correlation coefficient = 0.9994151
******************************************************************
Deviance= 815.9172
AIC= 1405.677
The distribution of the residuals is less heteroscedastic than before and the fit to the data better. The full summary of the model is:
******************************************************************
Family: c("NO", "Normal")
Call: gamlss(formula = SR ~ 1 + Sex + random(Text) + random(Language) + random(Speaker), sigma.formula = ~1 + Sex + random(Text) + random(Language) + random(Speaker), family = NO(mu.link = "identity"),
data = d, control = gamlss.control(n.cyc = 800, trace = FALSE), i.control = glim.control(bf.cyc = 800))
Fitting method: RS()
------------------------------------------------------------------
Mu link function: identity
Mu Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 6.461395 0.007705 838.58 <2e-16 ***
SexM 0.336314 0.011632 28.91 <2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
------------------------------------------------------------------
Sigma link function: log
Sigma Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.29481 0.02093 -61.853 < 2e-16 ***
SexM 0.11575 0.02957 3.915 9.34e-05 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
------------------------------------------------------------------
NOTE: Additive smoothing terms exist in the formulas:
i) Std. Error for smoothers are for the linear effect only.
ii) Std. Error for the linear terms maybe are not accurate.
------------------------------------------------------------------
No. of observations in the fit: 2288
Degrees of Freedom for the fit: 294.8797
Residual Deg. of Freedom: 1993.12
at cycle: 53
Global Deviance: 815.9172
AIC: 1405.677
SBC: 3096.939
******************************************************************
Text
Random effects fit using the gamlss function random()
Degrees of Freedom for the fit : 14.39185
Random effect parameter sigma_b: 0.109103
Smoothing parameter lambda : 84.5402
Language
Random effects fit using the gamlss function random()
Degrees of Freedom for the fit : 16.98344
Random effect parameter sigma_b: 0.870621
Smoothing parameter lambda : 1.32916
Speaker
Random effects fit using the gamlss function random()
Degrees of Freedom for the fit : 165.4211
Random effect parameter sigma_b: 0.569042
Smoothing parameter lambda : 3.32891
Text
Random effects fit using the gamlss function random()
Degrees of Freedom for the fit : 0.001525229
Random effect parameter sigma_b: 0.000518276
Smoothing parameter lambda : 3475490
Language
Random effects fit using the gamlss function random()
Degrees of Freedom for the fit : 14.62302
Random effect parameter sigma_b: 0.153321
Smoothing parameter lambda : 39.9687
Speaker
Random effects fit using the gamlss function random()
Degrees of Freedom for the fit : 79.45877
Random effect parameter sigma_b: 0.181484
Smoothing parameter lambda : 29.3639
******************************************************************
Summary of the Quantile Residuals
mean = 8.566117e-05
variance = 1.000437
coef. of skewness = 0.09721211
coef. of kurtosis = 3.684191
Filliben correlation coefficient = 0.9978017
******************************************************************
Deviance= 9322.636
AIC= 9721.821
******************************************************************
Summary of the Quantile Residuals
mean = 0.00210489
variance = 1.000442
coef. of skewness = 0.0313958
coef. of kurtosis = 2.864031
Filliben correlation coefficient = 0.9993924
******************************************************************
Deviance= 8961.553
AIC= 9554.279
Again, this is a better fit to the data. The full summary of the model is:
******************************************************************
Family: c("NO", "Normal")
Call: gamlss(formula = IR ~ 1 + Sex + random(Text) + random(Language) + random(Speaker), sigma.formula = ~1 + Sex + random(Text) + random(Language) + random(Speaker), family = NO(mu.link = "identity"),
data = d, control = gamlss.control(n.cyc = 800, trace = FALSE), i.control = glim.control(bf.cyc = 800))
Fitting method: RS()
------------------------------------------------------------------
Mu link function: identity
Mu Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 38.13387 0.04563 835.69 <2e-16 ***
SexM 2.02127 0.06926 29.18 <2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
------------------------------------------------------------------
Sigma link function: log
Sigma Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.48619 0.02093 23.225 < 2e-16 ***
SexM 0.11870 0.02957 4.015 6.17e-05 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
------------------------------------------------------------------
NOTE: Additive smoothing terms exist in the formulas:
i) Std. Error for smoothers are for the linear effect only.
ii) Std. Error for the linear terms maybe are not accurate.
------------------------------------------------------------------
No. of observations in the fit: 2288
Degrees of Freedom for the fit: 296.3631
Residual Deg. of Freedom: 1991.637
at cycle: 6
Global Deviance: 8961.553
AIC: 9554.279
SBC: 11254.05
******************************************************************
Text
Random effects fit using the gamlss function random()
Degrees of Freedom for the fit : 14.41202
Random effect parameter sigma_b: 0.660353
Smoothing parameter lambda : 2.30777
Language
Random effects fit using the gamlss function random()
Degrees of Freedom for the fit : 16.95235
Random effect parameter sigma_b: 3.10086
Smoothing parameter lambda : 0.104779
Speaker
Random effects fit using the gamlss function random()
Degrees of Freedom for the fit : 165.2917
Random effect parameter sigma_b: 3.39374
Smoothing parameter lambda : 0.0935855
Text
Random effects fit using the gamlss function random()
Degrees of Freedom for the fit : 0.0004688671
Random effect parameter sigma_b: 0.000282873
Smoothing parameter lambda : 11665300
Language
Random effects fit using the gamlss function random()
Degrees of Freedom for the fit : 14.73161
Random effect parameter sigma_b: 0.157608
Smoothing parameter lambda : 37.8204
Speaker
Random effects fit using the gamlss function random()
Degrees of Freedom for the fit : 80.97494
Random effect parameter sigma_b: 0.184921
Smoothing parameter lambda : 28.2979
Let’s model SR with ID as an additional predictor (fixed effect). N.B. In this case, we must drop Language as a random effect, since each language has, by definition, only one value of ID.
******************************************************************
Family: c("NO", "Normal")
Call: gamlss(formula = SR ~ 1 + ID + Sex + random(Text) + random(Speaker), sigma.formula = ~1 + ID + Sex + random(Text) + random(Speaker), family = NO(mu.link = "identity"), data = d, control = gamlss.control(n.cyc = 800,
trace = FALSE), i.control = glim.control(bf.cyc = 800))
Fitting method: RS()
------------------------------------------------------------------
Mu link function: identity
Mu Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 11.804255 0.037297 316.50 <2e-16 ***
ID -0.888578 0.005922 -150.06 <2e-16 ***
SexM 0.330217 0.011535 28.63 <2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
------------------------------------------------------------------
Sigma link function: log
Sigma Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.79327 0.10405 -7.624 3.77e-14 ***
ID -0.08460 0.01694 -4.993 6.45e-07 ***
SexM 0.11447 0.02957 3.872 0.000112 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
------------------------------------------------------------------
NOTE: Additive smoothing terms exist in the formulas:
i) Std. Error for smoothers are for the linear effect only.
ii) Std. Error for the linear terms maybe are not accurate.
------------------------------------------------------------------
No. of observations in the fit: 2288
Degrees of Freedom for the fit: 292.3552
Residual Deg. of Freedom: 1995.645
at cycle: 8
Global Deviance: 790.205
AIC: 1374.915
SBC: 3051.699
******************************************************************
******************************************************************
Summary of the Quantile Residuals
mean = 0.00189978
variance = 1.000433
coef. of skewness = 0.01948555
coef. of kurtosis = 2.773947
Filliben correlation coefficient = 0.9993363
******************************************************************
Deviance= 790.205
AIC= 1374.915
Adding ID as a predictor improves the fits (as judged by AIC). There is a negative estimate for ID, but significance is difficult to assess with GAMLSS model involving smoothing functions. However, also using a simple lmer model we have a significant effect of ID:
Type III Analysis of Variance Table with Satterthwaite's method
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
ID 16.8113 16.8113 1 166.64 158.2216 < 2.2e-16 ***
Sex 0.8857 0.8857 1 166.73 8.3361 0.004402 **
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Text
Random effects fit using the gamlss function random()
Degrees of Freedom for the fit : 14.41712
Random effect parameter sigma_b: 0.110654
Smoothing parameter lambda : 82.1886
Speaker
Random effects fit using the gamlss function random()
Degrees of Freedom for the fit : 167.487
Random effect parameter sigma_b: 0.757681
Smoothing parameter lambda : 1.8795
Text
Random effects fit using the gamlss function random()
Degrees of Freedom for the fit : 0.03164023
Random effect parameter sigma_b: 0.0023907
Smoothing parameter lambda : 155417
Speaker
Random effects fit using the gamlss function random()
Degrees of Freedom for the fit : 104.4195
Random effect parameter sigma_b: 0.242521
Smoothing parameter lambda : 15.8244
Linear mixed model fit by REML. t-tests use Satterthwaite's method ['lmerModLmerTest']
Formula: SR ~ Age * Sex + (1 | Text) + (1 | Language)
Data: info.rate.data
REML criterion at convergence: 3740.9
Scaled residuals:
Min 1Q Median 3Q Max
-3.4947 -0.6305 0.0105 0.5978 3.6970
Random effects:
Groups Name Variance Std.Dev.
Text (Intercept) 0.01299 0.1140
Language (Intercept) 1.01328 1.0066
Residual 0.36283 0.6024
Number of obs: 1979, groups: Text, 15; Language, 14
Fixed effects:
Estimate Std. Error df t value Pr(>|t|)
(Intercept) 6.701e+00 2.799e-01 1.516e+01 23.935 1.85e-13 ***
Age -7.195e-03 2.266e-03 1.950e+03 -3.175 0.00152 **
SexM 2.167e-01 9.216e-02 1.948e+03 2.351 0.01880 *
Age:SexM 2.446e-03 2.881e-03 1.949e+03 0.849 0.39597
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Correlation of Fixed Effects:
(Intr) Age SexM
Age -0.246
SexM -0.183 0.682
Age:SexM 0.176 -0.716 -0.956
Linear mixed model fit by REML. t-tests use Satterthwaite's method ['lmerModLmerTest']
Formula: IR ~ Age * Sex + (1 | Text) + (1 | Language)
Data: info.rate.data
REML criterion at convergence: 10789.9
Scaled residuals:
Min 1Q Median 3Q Max
-3.2281 -0.6258 0.0057 0.5889 4.6591
Random effects:
Groups Name Variance Std.Dev.
Text (Intercept) 0.4338 0.6587
Language (Intercept) 10.6184 3.2586
Residual 12.9843 3.6034
Number of obs: 1979, groups: Text, 15; Language, 14
Fixed effects:
Estimate Std. Error df t value Pr(>|t|)
(Intercept) 39.52601 0.98522 20.85678 40.119 < 2e-16 ***
Age -0.04694 0.01355 1954.05630 -3.465 0.000542 ***
SexM 1.34584 0.55127 1949.50697 2.441 0.014720 *
Age:SexM 0.01441 0.01723 1949.67588 0.836 0.403070
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Correlation of Fixed Effects:
(Intr) Age SexM
Age -0.418
SexM -0.310 0.683
Age:SexM 0.300 -0.716 -0.956
So, it seems * Age* and Sex are both worth including in our models (even if we have to discard quite a bit of data because of missing Age info). (In fact, the effect of Age seems more significant than that of Sex.)
In the following, we investigate if Age does matter when using GAMLSS modelling…
Bescause there is missing data fro Age, and because the GAMLSS models require no missing data, we will fit the models with Age (and its interaction with Sex) on the subset of the data that contains only those speakers with Age info. To make comparability possible, we also fit the same models but without Age on the exact same subset of the data.
******************************************************************
Summary of the Quantile Residuals
mean = -8.252317e-05
variance = 1.000506
coef. of skewness = 0.04585
coef. of kurtosis = 3.821224
Filliben correlation coefficient = 0.9973638
******************************************************************
The model including Age * Sex is:
******************************************************************
Family: c("NO", "Normal")
Call: gamlss(formula = SR ~ 1 + Sex * Age + random(Text) + random(Language) + random(Speaker), family = NO(mu.link = "identity"), data = info.rate.data.for.age, control = gamlss.control(n.cyc = 800,
trace = FALSE), i.control = glim.control(bf.cyc = 800))
Fitting method: RS()
------------------------------------------------------------------
Mu link function: identity
Mu Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 6.650290 0.034398 193.336 < 2e-16 ***
SexM 0.183376 0.046941 3.907 9.71e-05 ***
Age -0.004207 0.001077 -3.907 9.68e-05 ***
SexM:Age 0.003449 0.001464 2.356 0.0186 *
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
------------------------------------------------------------------
Sigma link function: log
Sigma Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.1605 0.0159 -73.01 <2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
------------------------------------------------------------------
NOTE: Additive smoothing terms exist in the formulas:
i) Std. Error for smoothers are for the linear effect only.
ii) Std. Error for the linear terms maybe are not accurate.
------------------------------------------------------------------
No. of observations in the fit: 1979
Degrees of Freedom for the fit: 161.6772
Residual Deg. of Freedom: 1817.323
at cycle: 35
Global Deviance: 1022.798
AIC: 1346.152
SBC: 2249.984
******************************************************************
The compared models are:
| Model | Deviance | AIC |
|---|---|---|
| Age * Sex | 1022.8 | 1346.2 |
| Age + Sex | 1022.8 | 1344.2 |
| Sex | 1022.8 | 1342.2 |
So, even if Age has a significant (negative) effect and interaction with Sex (positive for males), adding it does not seem to be warranted here…
******************************************************************
Summary of the Quantile Residuals
mean = 0.002353286
variance = 1.000499
coef. of skewness = 0.02299311
coef. of kurtosis = 2.90592
Filliben correlation coefficient = 0.9993823
******************************************************************
The model including Age * Sex is:
******************************************************************
Family: c("NO", "Normal")
Call: gamlss(formula = SR ~ 1 + Sex * Age + random(Text) + random(Language) + random(Speaker), sigma.formula = ~1 + Sex * Age + random(Text) + random(Language) + random(Speaker),
family = NO(mu.link = "identity"), data = info.rate.data.for.age, control = gamlss.control(n.cyc = 800, trace = FALSE), i.control = glim.control(bf.cyc = 800))
Fitting method: RS()
------------------------------------------------------------------
Mu link function: identity
Mu Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 6.5022917 0.0274165 237.167 < 2e-16 ***
SexM 0.0884641 0.0407203 2.172 0.030 *
Age 0.0007252 0.0008540 0.849 0.396
SexM:Age 0.0061826 0.0012966 4.768 2.01e-06 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
------------------------------------------------------------------
Sigma link function: log
Sigma Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.3105055 0.0798331 -16.416 <2e-16 ***
SexM -0.0396906 0.1079403 -0.368 0.713
Age 0.0004675 0.0025036 0.187 0.852
SexM:Age 0.0051227 0.0033693 1.520 0.129
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
------------------------------------------------------------------
NOTE: Additive smoothing terms exist in the formulas:
i) Std. Error for smoothers are for the linear effect only.
ii) Std. Error for the linear terms maybe are not accurate.
------------------------------------------------------------------
No. of observations in the fit: 1979
Degrees of Freedom for the fit: 243.8534
Residual Deg. of Freedom: 1735.147
at cycle: 10
Global Deviance: 721.7862
AIC: 1209.493
SBC: 2572.718
******************************************************************
The compared models are:
| Model | Deviance | AIC |
|---|---|---|
| Age * Sex | 721.8 | 1209.5 |
| Age + Sex | 721.5 | 1206.1 |
| Sex | 721.3 | 1203 |
So, even if Age has a significant (negative) effect (but no interaction with Sex), adding it does not seem to be warranted here either…
The distribution of the residuals is less heteroscedastic than before and the fit to the data better.
Thus, for SR, even if there is a hint that Age might affect it negatively (and there might also be an interaction with Sex with a positive effect for males), overall, the various fit indices do not warrant its inclusion in the GAMLSS models.
******************************************************************
Summary of the Quantile Residuals
mean = 1.384544e-05
variance = 1.000506
coef. of skewness = 0.09166289
coef. of kurtosis = 3.764844
Filliben correlation coefficient = 0.9974381
******************************************************************
The model including Age * Sex is:
******************************************************************
Family: c("NO", "Normal")
Call: gamlss(formula = IR ~ 1 + Sex * Age + random(Text) + random(Language) + random(Speaker), family = NO(mu.link = "identity"), data = info.rate.data.for.age, control = gamlss.control(n.cyc = 800,
trace = FALSE), i.control = glim.control(bf.cyc = 800))
Fitting method: RS()
------------------------------------------------------------------
Mu link function: identity
Mu Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 39.839840 0.201527 197.689 < 2e-16 ***
SexM 0.350589 0.275016 1.275 0.203
Age -0.058929 0.006309 -9.341 < 2e-16 ***
SexM:Age 0.046375 0.008576 5.408 7.22e-08 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
------------------------------------------------------------------
Sigma link function: log
Sigma Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.6074 0.0159 38.21 <2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
------------------------------------------------------------------
NOTE: Additive smoothing terms exist in the formulas:
i) Std. Error for smoothers are for the linear effect only.
ii) Std. Error for the linear terms maybe are not accurate.
------------------------------------------------------------------
No. of observations in the fit: 1979
Degrees of Freedom for the fit: 161.8269
Residual Deg. of Freedom: 1817.173
at cycle: 2
Global Deviance: 8020.303
AIC: 8343.956
SBC: 9248.625
******************************************************************
The compared models are:
| Model | Deviance | AIC |
|---|---|---|
| Age * Sex | 8020.3 | 8344 |
| Age + Sex | 8020.3 | 8342 |
| Sex | 8020.3 | 8340 |
So, even if Age has a significant (negative) effect and interaction with Sex (positive for males) – interestingly, in this case the main effect of Sex disappears –, adding it does not seem to be warranted…
******************************************************************
Summary of the Quantile Residuals
mean = 0.002387767
variance = 1.00051
coef. of skewness = 0.02686616
coef. of kurtosis = 2.916301
Filliben correlation coefficient = 0.9993496
******************************************************************
The model including Age * Sex is:
******************************************************************
Family: c("NO", "Normal")
Call: gamlss(formula = IR ~ 1 + Sex * Age + random(Text) + random(Language) + random(Speaker), sigma.formula = ~1 + Sex * Age + random(Text) + random(Language) + random(Speaker),
family = NO(mu.link = "identity"), data = info.rate.data.for.age, control = gamlss.control(n.cyc = 800, trace = FALSE), i.control = glim.control(bf.cyc = 800))
Fitting method: RS()
------------------------------------------------------------------
Mu link function: identity
Mu Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 39.742235 0.159812 248.681 < 2e-16 ***
SexM 0.598655 0.236910 2.527 0.0116 *
Age -0.055990 0.004891 -11.449 < 2e-16 ***
SexM:Age 0.038192 0.007416 5.150 2.9e-07 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
------------------------------------------------------------------
Sigma link function: log
Sigma Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.5044688 0.0797812 6.323 3.25e-10 ***
SexM -0.0500481 0.1079793 -0.463 0.643
Age -0.0008933 0.0025018 -0.357 0.721
SexM:Age 0.0054948 0.0033706 1.630 0.103
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
------------------------------------------------------------------
NOTE: Additive smoothing terms exist in the formulas:
i) Std. Error for smoothers are for the linear effect only.
ii) Std. Error for the linear terms maybe are not accurate.
------------------------------------------------------------------
No. of observations in the fit: 1979
Degrees of Freedom for the fit: 244.3666
Residual Deg. of Freedom: 1734.633
at cycle: 10
Global Deviance: 7726.654
AIC: 8215.387
SBC: 9581.481
******************************************************************
The compared models are:
| Model | Deviance | AIC |
|---|---|---|
| Age * Sex | 7726.7 | 8215.4 |
| Age + Sex | 7726.4 | 8212.1 |
| Sex | 7725.8 | 8208.7 |
So, even if Age has a significant (negative) effect and interaction with Sex (positive for males) – interestingly, in this case the main effect of Sex disappears –, adding it does not seem to be warranted…
The distribution of the residuals is less heteroscedastic than before and the fit to the data better.
Thus, while for IR the hint that Age has a negative main effect and interacts with Sex (with a positive effect for males, containing the whole effect of Sex) is much stronger, the various fit indices do not warrant its inclusion in the GAMLSS models.
******************************************************************
Summary of the Quantile Residuals
mean = 0.001744233
variance = 1.000502
coef. of skewness = 0.008063235
coef. of kurtosis = 2.834363
Filliben correlation coefficient = 0.9993747
******************************************************************
The model including Age * Sex is:
******************************************************************
Family: c("NO", "Normal")
Call: gamlss(formula = SR ~ 1 + ID + Sex * Age + random(Text) + random(Speaker), sigma.formula = ~1 + ID + Sex * Age + random(Text) + random(Speaker), family = NO(mu.link = "identity"),
data = info.rate.data.for.age, control = gamlss.control(n.cyc = 800, trace = FALSE), i.control = glim.control(bf.cyc = 800))
Fitting method: RS()
------------------------------------------------------------------
Mu link function: identity
Mu Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 12.6896560 0.0553430 229.291 < 2e-16 ***
ID -0.9813601 0.0071683 -136.904 < 2e-16 ***
SexM -0.0205028 0.0408781 -0.502 0.616
Age -0.0109959 0.0008679 -12.670 < 2e-16 ***
SexM:Age 0.0100110 0.0013074 7.657 3.12e-14 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
------------------------------------------------------------------
Sigma link function: log
Sigma Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.743754 0.145088 -5.126 3.29e-07 ***
ID -0.086079 0.019128 -4.500 7.24e-06 ***
SexM -0.092987 0.108736 -0.855 0.3926
Age -0.001499 0.002518 -0.595 0.5518
SexM:Age 0.006919 0.003396 2.037 0.0418 *
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
------------------------------------------------------------------
NOTE: Additive smoothing terms exist in the formulas:
i) Std. Error for smoothers are for the linear effect only.
ii) Std. Error for the linear terms maybe are not accurate.
------------------------------------------------------------------
No. of observations in the fit: 1979
Degrees of Freedom for the fit: 236.2345
Residual Deg. of Freedom: 1742.766
at cycle: 8
Global Deviance: 705.2124
AIC: 1177.681
SBC: 2498.314
******************************************************************
The compared models are:
| Model | Deviance | AIC |
|---|---|---|
| Age * Sex | 705.2 | 1177.7 |
| Age + Sex | 704.5 | 1174.4 |
| Sex | 704.6 | 1170.8 |
Clearly, adding Age is not warranted here…
As above, we also looked a the simple lmer model:
The compared models are:
| Model | AIC |
|---|---|
| Age * Sex | 1818.2 |
| Age + Sex | 1809.9 |
| Sex | 1800.4 |
Interestingly, here the best model (as suggested by AIC) is the one with both Age and Sex (but no interaaction), but, as can be seen below, actually the effect of Age is not significant (p = 0.402), suggesting that, in fact, the best model is still the one not including Age:
Linear mixed model fit by REML. t-tests use Satterthwaite's method ['lmerModLmerTest']
Formula: SR ~ 1 + ID + Sex + Age + (1 | Text) + (1 | Speaker)
Data: info.rate.data.for.age
REML criterion at convergence: 1795.9
Scaled residuals:
Min 1Q Median 3Q Max
-3.8108 -0.6235 0.0128 0.5918 5.1698
Random effects:
Groups Name Variance Std.Dev.
Speaker (Intercept) 0.52573 0.7251
Text (Intercept) 0.01494 0.1222
Residual 0.10575 0.3252
Number of obs: 1979, groups: Speaker, 132; Text, 15
Fixed effects:
Estimate Std. Error df t value Pr(>|t|)
(Intercept) 12.526462 0.537328 128.871889 23.313 <2e-16 ***
ID -0.981475 0.076625 128.000256 -12.809 <2e-16 ***
SexM 0.286423 0.127080 128.001446 2.254 0.0259 *
Age -0.005613 0.006674 128.000190 -0.841 0.4019
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Correlation of Fixed Effects:
(Intr) ID SexM
ID -0.910
SexM -0.114 -0.005
Age -0.521 0.167 -0.005
Within the limits of our reduced dataset (containing only 132 speakers with Age info), we found the following:
When modelling SR and IR with GAMLSS, while there are hints that Age has, for both, overall:
it does not seem warranted to include it in these models.
When modelling the relationship between SR and ID, this negative relationship:
but, alas, the inclusion of Age is not warranted in the GAMLSS model, nor (really) in the simpler LMER model.
Thus, while Age seems to negatively influence (in a sex-dependent manner) both SR and IR, as well as strengthen the negative relationship between them, its effects are far from clear in the current dataset.
In what follows, mixing probabilities are independent from factors such as Sex.
Between 1 and 5 Gaussian distributions:
1 component
Mixing Family: "NO"
Fitting method: EM algorithm
Call: gamlssMX(formula = SR ~ 1, family = NO, K = 1, data = d, plot = FALSE)
Mu Coefficients for model: 1
(Intercept)
6.631
Sigma Coefficients for model: 1
(Intercept)
0.1378
Estimated probabilities: 1
Degrees of Freedom for the fit: 2 Residual Deg. of Freedom 2286
Global Deviance: 7123.61
AIC: 7127.61
SBC: 7139.08
2 components
Mixing Family: c("NO", "NO")
Fitting method: EM algorithm
Call: gamlssMX(formula = SR ~ 1, family = NO, K = 2, data = d, plot = FALSE)
Mu Coefficients for model: 1
(Intercept)
7.149
Sigma Coefficients for model: 1
(Intercept)
-0.1855
Mu Coefficients for model: 2
(Intercept)
5.259
Sigma Coefficients for model: 2
(Intercept)
-0.4788
Estimated probabilities: 0.7260515 0.2739485
Degrees of Freedom for the fit: 5 Residual Deg. of Freedom 2283
Global Deviance: 7001.22
AIC: 7011.22
SBC: 7039.9
3 components
Mixing Family: c("NO", "NO", "NO")
Fitting method: EM algorithm
Call: gamlssMX(formula = SR ~ 1, family = NO, K = 3, data = d, plot = FALSE)
Mu Coefficients for model: 1
(Intercept)
7.106
Sigma Coefficients for model: 1
(Intercept)
-0.2769
Mu Coefficients for model: 2
(Intercept)
7.207
Sigma Coefficients for model: 2
(Intercept)
-0.133
Mu Coefficients for model: 3
(Intercept)
5.27
Sigma Coefficients for model: 3
(Intercept)
-0.4741
Estimated probabilities: 0.3374419 0.3826253 0.2799327
Degrees of Freedom for the fit: 8 Residual Deg. of Freedom 2280
Global Deviance: 7000.8
AIC: 7016.8
SBC: 7062.69
4 components
Mixing Family: c("NO", "NO", "NO", "NO")
Fitting method: EM algorithm
Call: gamlssMX(formula = SR ~ 1, family = NO, K = 4, data = d, plot = FALSE)
Mu Coefficients for model: 1
(Intercept)
7.784
Sigma Coefficients for model: 1
(Intercept)
-0.3468
Mu Coefficients for model: 2
(Intercept)
5.811
Sigma Coefficients for model: 2
(Intercept)
-0.2546
Mu Coefficients for model: 3
(Intercept)
7.124
Sigma Coefficients for model: 3
(Intercept)
-0.7752
Mu Coefficients for model: 4
(Intercept)
5.442
Sigma Coefficients for model: 4
(Intercept)
-0.3356
Estimated probabilities: 0.2584335 0.2204151 0.2985522 0.2225992
Degrees of Freedom for the fit: 11 Residual Deg. of Freedom 2277
Global Deviance: 6990.94
AIC: 7012.94
SBC: 7076.03
5 components
Mixing Family: c("NO", "NO", "NO", "NO", "NO")
Fitting method: EM algorithm
Call: gamlssMX(formula = SR ~ 1, family = NO, K = 5, data = d, plot = FALSE)
Mu Coefficients for model: 1
(Intercept)
6.396
Sigma Coefficients for model: 1
(Intercept)
-0.3948
Mu Coefficients for model: 2
(Intercept)
6.504
Sigma Coefficients for model: 2
(Intercept)
-0.3149
Mu Coefficients for model: 3
(Intercept)
5.132
Sigma Coefficients for model: 3
(Intercept)
-0.5414
Mu Coefficients for model: 4
(Intercept)
7.267
Sigma Coefficients for model: 4
(Intercept)
-1.037
Mu Coefficients for model: 5
(Intercept)
7.959
Sigma Coefficients for model: 5
(Intercept)
-0.4526
Estimated probabilities: 0.2027298 0.180915 0.2175461 0.191181 0.2076281
Degrees of Freedom for the fit: 14 Residual Deg. of Freedom 2274
Global Deviance: 6981.9
AIC: 7009.9
SBC: 7090.2
Comparing AIC
df AIC
mix.SR.NO.5 14 7009.904
mix.SR.NO.2 5 7011.223
mix.SR.NO.4 11 7012.941
mix.SR.NO.3 8 7016.804
mix.SR.NO.1 2 7127.609
Showing the distributions
Mixture of Gaussians for SR.
Between 1 and 5 Gaussian distributions:
1 component
Mixing Family: "NO"
Fitting method: EM algorithm
Call: gamlssMX(formula = IR ~ 1, family = NO, K = 1, data = d, plot = FALSE)
Mu Coefficients for model: 1
(Intercept)
39.15
Sigma Coefficients for model: 1
(Intercept)
1.629
Estimated probabilities: 1
Degrees of Freedom for the fit: 2 Residual Deg. of Freedom 2286
Global Deviance: 13945.2
AIC: 13949.2
SBC: 13960.7
2 components
Mixing Family: c("NO", "NO")
Fitting method: EM algorithm
Call: gamlssMX(formula = IR ~ 1, family = NO, K = 2, data = d, plot = FALSE)
Mu Coefficients for model: 1
(Intercept)
41.27
Sigma Coefficients for model: 1
(Intercept)
1.863
Mu Coefficients for model: 2
(Intercept)
38.38
Sigma Coefficients for model: 2
(Intercept)
1.447
Estimated probabilities: 0.266011 0.733989
Degrees of Freedom for the fit: 5 Residual Deg. of Freedom 2283
Global Deviance: 13894.8
AIC: 13904.8
SBC: 13933.5
3 components
Mixing Family: c("NO", "NO", "NO")
Fitting method: EM algorithm
Call: gamlssMX(formula = IR ~ 1, family = NO, K = 3, data = d, plot = FALSE)
Mu Coefficients for model: 1
(Intercept)
40
Sigma Coefficients for model: 1
(Intercept)
0.8341
Mu Coefficients for model: 2
(Intercept)
42.1
Sigma Coefficients for model: 2
(Intercept)
1.725
Mu Coefficients for model: 3
(Intercept)
36.09
Sigma Coefficients for model: 3
(Intercept)
1.395
Estimated probabilities: 0.2444004 0.350947 0.4046526
Degrees of Freedom for the fit: 8 Residual Deg. of Freedom 2280
Global Deviance: 13875.8
AIC: 13891.8
SBC: 13937.6
4 components
Mixing Family: c("NO", "NO", "NO", "NO")
Fitting method: EM algorithm
Call: gamlssMX(formula = IR ~ 1, family = NO, K = 4, data = d, plot = FALSE)
Mu Coefficients for model: 1
(Intercept)
41
Sigma Coefficients for model: 1
(Intercept)
1.127
Mu Coefficients for model: 2
(Intercept)
39.56
Sigma Coefficients for model: 2
(Intercept)
0.08543
Mu Coefficients for model: 3
(Intercept)
36.14
Sigma Coefficients for model: 3
(Intercept)
1.38
Mu Coefficients for model: 4
(Intercept)
42.51
Sigma Coefficients for model: 4
(Intercept)
1.774
Estimated probabilities: 0.222081 0.07674445 0.4389749 0.2621997
Degrees of Freedom for the fit: 11 Residual Deg. of Freedom 2277
Global Deviance: 13868.4
AIC: 13890.4
SBC: 13953.5
5 components
Mixing Family: c("NO", "NO", "NO", "NO", "NO")
Fitting method: EM algorithm
Call: gamlssMX(formula = IR ~ 1, family = NO, K = 5, data = d, plot = FALSE)
Mu Coefficients for model: 1
(Intercept)
40.91
Sigma Coefficients for model: 1
(Intercept)
1.303
Mu Coefficients for model: 2
(Intercept)
34.42
Sigma Coefficients for model: 2
(Intercept)
1.252
Mu Coefficients for model: 3
(Intercept)
39.47
Sigma Coefficients for model: 3
(Intercept)
0.5634
Mu Coefficients for model: 4
(Intercept)
43.68
Sigma Coefficients for model: 4
(Intercept)
1.767
Mu Coefficients for model: 5
(Intercept)
39.64
Sigma Coefficients for model: 5
(Intercept)
1.289
Estimated probabilities: 0.211754 0.2798107 0.1459857 0.1801956 0.182254
Degrees of Freedom for the fit: 14 Residual Deg. of Freedom 2274
Global Deviance: 13868.7
AIC: 13896.7
SBC: 13977
Comparing AIC
df AIC
mix.IR.NO.4 11 13890.37
mix.IR.NO.3 8 13891.76
mix.IR.NO.5 14 13896.71
mix.IR.NO.2 5 13904.80
mix.IR.NO.1 2 13949.21
Showing the distributions
Mixture of Gaussians for IR.
We used three ways to estimate how unimodal a distribution is, as they tend to disagree and the problem of unimodality testing is far from settled (see Freeman & Dale, 2013):
diptest;For each such test, we performed four randomisation procedures to obtain an estimate of the “specialness” of the observed unimodality estimate; for each new permuted dataset, we recompute everything before estimating the unimodlaity of the permuted distribution:
The observed estimate (vertical blue solid line), the permuted distribution (gray histogram), and the “unimodality region” (shaded green rectangle) are shown below (for PM3, we also show the original estimate using the Speaker average SR as a vertical solid red line).
Permutation of the texts’ SRs (PM1).
Permutation of the languages’ ID (PM2).
Permutation of the speakers’ average SRs (PM3).
Permutation of the languages’ average SRs with speaker adjustement (PM4).
| Scenario | Measure | Test | Observed estimate (p-value) | % more unimodal permutations |
|---|---|---|---|---|
| PM1 | SR | Silverman | - (0.017) | 98.4% |
| PM1 | SR | Dip | 0.005 (0.984) * | 100% |
| PM1 | SR | BC | 0.19 () * | 100% |
| PM1 | IR | Silverman | - (0.85) * | 16.5% |
| PM1 | IR | Dip | 0.005 (0.992) * | 72.5% |
| PM1 | IR | BC | 0.167 () * | 100% |
| PM2 | SR | Silverman | - (0.017) | 98.3% |
| PM2 | SR | Dip | 0.005 (0.984) * | 100% |
| PM2 | SR | BC | 0.19 () * | 100% |
| PM2 | IR | Silverman | - (0.85) * | 2.3% |
| PM2 | IR | Dip | 0.005 (0.992) * | 27.5% |
| PM2 | IR | BC | 0.167 () * | 96.3% |
| PM3 | SR | Silverman | - (0.017) | 100% |
| PM3 | SR | Dip | 0.005 (0.984) * | 0% |
| PM3 | SR | BC | 0.19 () * | 100% |
| PM3 | IR | Silverman | - (0.85) * | 14.9% |
| PM3 | IR | Dip | 0.005 (0.992) * | 16.9% |
| PM3 | IR | BC | 0.167 () * | 99.9% |
| PM4 | SR | Silverman | - (0.017) | 7.5% |
| PM4 | SR | Dip | 0.005 (0.984) * | 3.2% |
| PM4 | SR | BC | 0.19 () * | 84.3% |
| PM4 | IR | Silverman | - (0.85) * | 0.7% |
| PM4 | IR | Dip | 0.005 (0.992) * | 14.6% |
| PM4 | IR | BC | 0.167 () * | 97.7% |
We compute various distances between languages (as implemented by function distance() in package philentropy) in what concerns the distribution of NS, SR and ID.
Comparing the distribution of pairwise distances between languages.
| m1 | m2 | d | mean1 | median1 | sd1 | mean2 | median2 | sd2 | p |
|---|---|---|---|---|---|---|---|---|---|
| IR | NS | Hellinger | 0.88 | 0.83 | 0.32 | 1.20 | 1.19 | 0.39 | 0.00 |
| IR | NS | Jensen-Shannon | 0.17 | 0.14 | 0.12 | 0.29 | 0.26 | 0.16 | 0.00 |
| IR | NS | Kolmogorov–Smirnov | 0.42 | 0.37 | 0.20 | 0.57 | 0.57 | 0.23 | 0.00 |
| IR | NS | Kullback-Leibler | 7.13 | 4.22 | 7.81 | 15.42 | 13.12 | 13.31 | 0.00 |
| IR | NS | Squared-Chi | 0.56 | 0.45 | 0.36 | 0.88 | 0.79 | 0.47 | 0.00 |
| IR | SR | Hellinger | 0.88 | 0.83 | 0.32 | 1.10 | 1.06 | 0.48 | 0.00 |
| IR | SR | Jensen-Shannon | 0.17 | 0.14 | 0.12 | 0.27 | 0.23 | 0.20 | 0.00 |
| IR | SR | Kolmogorov–Smirnov | 0.42 | 0.37 | 0.20 | 0.56 | 0.55 | 0.27 | 0.00 |
| IR | SR | Kullback-Leibler | 7.13 | 4.22 | 7.81 | 12.80 | 6.69 | 14.88 | 0.00 |
| IR | SR | Squared-Chi | 0.56 | 0.45 | 0.36 | 0.86 | 0.77 | 0.57 | 0.00 |
| NS | SR | Hellinger | 1.20 | 1.19 | 0.39 | 1.10 | 1.06 | 0.48 | 0.01 |
| NS | SR | Jensen-Shannon | 0.29 | 0.26 | 0.16 | 0.27 | 0.23 | 0.20 | 0.26 |
| NS | SR | Kolmogorov–Smirnov | 0.57 | 0.57 | 0.23 | 0.56 | 0.55 | 0.27 | 0.83 |
| NS | SR | Kullback-Leibler | 15.42 | 13.12 | 13.31 | 12.80 | 6.69 | 14.88 | 0.05 |
| NS | SR | Squared-Chi | 0.88 | 0.79 | 0.47 | 0.86 | 0.77 | 0.57 | 0.60 |
Campione, E., & Véronis, J. (1998). A multilingual prosodic database, Proc. of the 5th International Conference on Spoken Language Pro cessing (ICSLP’98), Sydney, Australia, 3163-3166.
Freeman, J. B., & Dale, R. (2013). Assessing bimodality to detect the presence of a dual cognitive process. Behavior research methods, 45(1), 83-97.
Hall, P., & York, M. (2001). On the calibration of Silverman’s test for multimodality. Statistica Sinica, 11, 515-536.
Hartigan, J. A., & Hartigan, P. M. (1985) The Dip Test of Unimodality. Annals of Statistics 13, 70–84.
Hartigan, P. M. (1985) Computation of the Dip Statistic to Test for Unimodality. Applied Statistics (JRSS C) 34, 320–325.
Le, V. B., Tran, D. D., Castelli, E., Besacier, L., & Serignat, J. F. (2004). Spoken and Written Language Resources for Vietnamese. In LREC. 4, pp. 599-602.
Lyding, V., Stemle, E., Borghetti, C., Brunello, M., Castagnoli, S., Dell’Orletta, F., Dittmann, H., Lenci, A., & Pirrelli, V. (2014). The PAISÀ Corpus of Italian Web Texts. In Proceedings of the 9th Web as Corpus Workshop (WaC-9). Association for Computational Linguistics, Gothenburg, Sweden, 36-43.
New B., Pallier C., Ferrand L., & Matos R. (2001). Une base de données lexicales du français contemporain sur internet: LEXIQUE 3.80, L’Année Psychologique, 101, 447-462. http://www.lexique.org.
Oh, Y. M. (2015). Linguistic complexity and information: quantitative approaches. PhD Thesis, Université de Lyon, France. Retrieved from http://www.afcp-parole.org/doc/theses/these_YMO15.pdf
Perea, M., Urkia, M., Davis, C. J., Agirre, A., Laseka, E., & Carreiras, M. (2006). E-Hitz: A word frequency list and a program for deriving psycholinguistic statistics in an agglutinative language (Basque). Behavior Research Methods, 38(4), 610-615.
Sharoff, S. (2006). Creating general-purpose corpora using automated search engine queries. In Baroni, M. and Bernardini, S. (Eds.) WaCky! Working papers on the web as corpus, Gedit, Bologna, http://corpus.leeds.ac.uk/queryzh.html.
Silverman, B.W. (1981). Using Kernel Density Estimates to investigate Multimodality. Journal of the Royal Statistical Society, Series B, 43, 97-99.
Váradi, T. (2002). The Hungarian National Corpus. In LREC.
Zséder, A., Recski, G., Varga, D., & Kornai, A. (2012). Rapid creation of large-scale corpora and frequency dictionaries. In Proceedings to LREC 2012.
R session infoThis document was compiled on:
R version 3.4.1 (2017-06-30)
**Platform:** x86_64-pc-linux-gnu (64-bit)
locale: LC_CTYPE=en_US.UTF-8, LC_NUMERIC=C, LC_TIME=en_US.UTF-8, LC_COLLATE=en_US.UTF-8, LC_MONETARY=en_US.UTF-8, LC_MESSAGES=en_US.UTF-8, LC_PAPER=en_US.UTF-8, LC_NAME=C, LC_ADDRESS=C, LC_TELEPHONE=C, LC_MEASUREMENT=en_US.UTF-8 and LC_IDENTIFICATION=C
attached base packages: compiler, grid, parallel, splines, stats, graphics, grDevices, utils, datasets, methods and base
other attached packages: broman(v.0.68-2), philentropy(v.0.3.0), diptest(v.0.75-7), pander(v.0.6.3), moments(v.0.14), sjPlot(v.2.6.2), sjstats(v.0.17.3), gamlss.mx(v.4.3-5), nnet(v.7.3-12), gamlss(v.5.1-2), nlme(v.3.1-137), gamlss.dist(v.5.1-1), MASS(v.7.3-51.1), gamlss.data(v.5.1-0), lmerTest(v.3.1-0), lme4(v.1.1-20), Matrix(v.1.2-15), plyr(v.1.8.4), reshape2(v.1.4.3), ggplot2(v.3.1.0) and RhpcBLASctl(v.0.18-205)
loaded via a namespace (and not attached): RColorBrewer(v.1.1-2), numDeriv(v.2016.8-1), tools(v.3.4.1), TMB(v.1.7.15), backports(v.1.1.3), R6(v.2.4.0), sjlabelled(v.1.0.16), lazyeval(v.0.2.1), colorspace(v.1.4-0), withr(v.2.1.2), tidyselect(v.0.2.5), mnormt(v.1.5-5), emmeans(v.1.3.2), sandwich(v.2.5-0), labeling(v.0.3), scales(v.1.0.0), mvtnorm(v.1.0-8), psych(v.1.8.12), ggridges(v.0.5.1), stringr(v.1.4.0), digest(v.0.6.18), foreign(v.0.8-71), minqa(v.1.2.4), rmarkdown(v.1.11), stringdist(v.0.9.5.1), pkgconfig(v.2.0.2), htmltools(v.0.3.6), highr(v.0.7), pwr(v.1.2-2), rlang(v.0.3.1), generics(v.0.0.2), zoo(v.1.8-4), dplyr(v.0.8.0.1), magrittr(v.1.5), modeltools(v.0.2-22), bayesplot(v.1.6.0), Rcpp(v.1.0.0), munsell(v.0.5.0), prediction(v.0.3.6.2), stringi(v.1.3.1), multcomp(v.1.4-8), yaml(v.2.2.0), snakecase(v.0.9.2), sjmisc(v.2.7.7), forcats(v.0.4.0), crayon(v.1.3.4), lattice(v.0.20-38), ggeffects(v.0.8.0), haven(v.2.1.0), hms(v.0.4.2), knitr(v.1.21), pillar(v.1.3.1), estimability(v.1.3), codetools(v.0.2-16), stats4(v.3.4.1), glue(v.1.3.0), evaluate(v.0.13), data.table(v.1.12.0), modelr(v.0.1.4), nloptr(v.1.2.1), gtable(v.0.2.0), purrr(v.0.3.0), tidyr(v.0.8.2), assertthat(v.0.2.0), xfun(v.0.5), coin(v.1.2-2), xtable(v.1.8-3), broom(v.0.5.1), coda(v.0.19-2), survival(v.2.43-3), tibble(v.2.0.1), glmmTMB(v.0.2.3) and TH.data(v.1.0-10)
Here we generate the figures used in the main paper (saved to the ./figures folder as 600 DPI TIFF files Figure-*.tiff).